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We study the effects of intersite RKKY-like interactions in the mixed valence 
regime of a two-band extended Hubbard modeh This model is known to display a 
metallic non-Fermi liquid state with spin-charge separation in the standard D = oo 
limit (where interactions lead to purely on-site fluctuations only) . Using an extended 
D = oo approach, we find that the spin-charge separation survives the quantum fluc- 
tuations associated with arbitrary intersite density-density interactions and a finite 
range of intersite spin-exchange interactions. We determine the qualitative behavior 
of the spin, charge, and single-particle correlation functions. The implications of 
these results for the solution in finite spatial dimensions are discussed. 
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Whether or not spin-charge separation exists in two or three dimensions is an important 
outstanding question in the area of strongly correlated electron systems. Spin-charge sepa- 
ration occurs in the Luttinger liquid, which describes interacting fermion systems in ID [|I|. 
The Luttinger-liquid picture breaks down in D > 1 when interactions are treated perturba- 
tively 0. What happens when interactions are treated non-perturbatively remains an open 
question Recently, in the opposite D = oo limit a non- Fermi liquid state with spin- 
charge separation has been identified in the mixed-valence regime of the two-band extended 
Hubbard model Named the intermediate phase, this state has quasiparticle-like spin 

excitations and incoherent single particle and charge excitations. 

A natural and important question at this point is whether the spin-charge-separation 
phenomenon persists in finite dimensions. Here, we address this issue using the extended 
D = oo approach @J^. The key observation motivating this approach is that, the most 
important ingredients that exist in finite dimensions but are missing in the standard D = oo 
limit are the quantum fiuctuations associated with intersite interactions. By introducing an 
appropriate scaling of the intersite two-particle interactions in terms of 1/D, the extended 
D = oo approach treats the intersite fiuctuations on an equal footing with on-site interac- 
tions. We note that several alternative approaches ||ID| , pT]] have also been proposed to go 
beyond the standard D = oo dynamical mean field theory; they appear to be less practical 
at this stage. 

The two-band extended Hubbard model we study is defined by the following Hamiltonian, 

i <ij> 

hi = t^Ud^ + Urid^^nd^i + Y,ii4aCia + H.c.) + Vud^Uc^ + JxSd, ■ s^^ 

a 

^ij = tij (^Lcja + Vij : Ud^ : : Ud^ : +JijSd, ■ Sd^ (1) 

There are two kinds of electrons, the strongly correlated d— electrons and the band 
c— electrons. Both have spin 1/2, with n^-, ric-, 5*^- and Sc- representing the correspond- 
ing density and spin operators at the site i. The Hamiltonian is written as the sum of 
two parts, hi describes purely local couplings at the site i. e° and U are the energy level 
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and the Hubbard interaction respectively. For the purpose of studying the mixed-valence 
regime it suffices to take U to be infinite, t is the hybridization. Also included are two more 
local interaction terms allowed by symmetry (which are neglected in the standard Anderson 
lattice model): ^ is a charge-screening interaction and Jk is a spin-exchange interaction. 
hij describes three kinds of intersite couplings, including the single-particle hopping (tjj), 
the two-particle RKKY-like density- density interaction (vij) and spin-exchange interaction 

The extended large D limit is taken by scaling the two-particle nearest-neighbor inter- 
action terms f<jj> = vo/V2D and J<ij> = J0/V2D, along with the usual scaling for the 
single-particle hopping t<jj> = to/V^D, while keeping vq, Jq, and to fixed. This limit is well 
defined when the static (Hartree) terms are subtracted by introducing the normal ordering: 
: n := n— < n >. We will limit our discussion to the non-ordered states. In the D = 00 
limit, all the correlation functions of the lattice model can be calculated from an impurity 
coupled to self-consistent media. The effective impurity action can be derived using the 
"cavity method" Ij^,^. One first selects an arbitrary site on the lattice and expands the 
partition function in terms of all non-zero hio + /iQj, each of which is of order One 
then integrates out all the degrees of freedom at sites other than site 0, and takes the D = 00 
limit. The resulting effective impurity action is as follows, 

Sili = So- dTdr'l^cl^ir) Go\t - r ) Co,^(r') 

+ -ndo- (r) Xci - -r') ■■ rido ■ (r ) + Sd^ir) x~o(^ " ' Sdo{r')] (2) 

where 5*0 is the action associated with ho, and Gq^, XcOj x7o determined self- 
consistently by the single particle Green function, the two-particle density-density and spin- 
spin correlation functions, respectively. The impurity action can be equivalently written 
in terms of the following impurity Hamiltonian, 

Himp = Hisin + E^ndf) + Undg-^ndoi + Vudorico + JxSdo ■ 

+tY,{dl,Ck, + H.c.) + Fj2:nd,:{p, + pi,) + gJ2 4 " (0, + 0-.) (3) 

ka q q 
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Eq. describes an impurity coupled to three non-interacting bands. The impurity has 
an energy level = — fi, where /i is the chemical potential, and a Hubbard interaction 
U. The three bands are, respectively, fermionic, scalar-bosonic, and vector-bosonic. We 
stress that the different components of the vector bosons commute with each other ||12|| : 
l^q^ ^q'^] — ^tiu^qq' whcrc ^,i> = X, y , z. The dispersions of these bands are specified by 

Hkin = J2kaEkcl„Cka + Hq^qPlPq + llqWq4>\ ' <Pq- Wc haVC choSCU Efc, Wq, aud Wq SUCh 

that the couplings of the impurity to the three bands are independent of momentum. In 
the mixed-valence regime, the self-consistent density of states of the fermionic bath at the 
Fermi energy, po? turns out to be always finite. This is a crucial feature that allows for the 
following asymptotic analysis. Without loss of generality, we will in addition consider only 
the case of the Bethe lattice with infinite connectivity. Here, the self-consistency equations 
for the scalar- and vector- bosonic baths have the following forms: 

/Ef-j^-^o^X.O^J (4) 

where Xc and Xs are the local dynamical charge and spin susceptibilities, respectively. 

When intersite couplings (Jo and vq) are absent, F and g vanish in H^JJ^. The corre- 
sponding self-consistent problem has been solved asymptotically exactly When Jk is 
antiferromagnetic, there is a phase transition between the usual strong-coupling phase, a 
Fermi liquid, and the intermediate phase, a non-Fermi liquid with spin-charge separation. 
In the presence of intersite couplings, both F and g become non-zero. The impurity d or- 
bital is then coupled to the additional scalar-bosonic and vector-bosonic baths. Two effects 
of these additional baths have to be addressed. First, do they affect the coupling between 
the spin and charge degrees of freedom such that the spin-charge-separation phenomenon is 
destroyed? Second, if spin-charge separation persists, what are the forms of the spin, charge, 
and single-particle correlation functions? 

To address these questions, we need first to know the solution to the impurity Hamilto- 
nian, H^JJ^^ for generic forms of Wq and Wq {Ek is such that po is always finite). Anticipating 



the scale-invariant forms for the charge and spin susceptibihties, we consider the case when 
Wg and Wg are such that, 
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Here ac, a^, Kc, and Kg are certain fixed parameters. (The self-consistency, Eq. @, is not 
imposed at this point). 

When g = [|l^], the effect of the F— coupling can be determined using a renormalization 
group (RG) analysis based on a Coulomb-gas representation for the partition function Zj^p = 
Tr(e~ In a standard fashion, we expand Zimp in terms of the t— and Jk~ couplings 

(where is the transverse component of the Ji<-— coupling) leading to an infinite series. 
Each term in the series is associated with a particular history of the impurity configurations, 
|0 >d and |cr >d= dllO >d, along the imaginary time axis [0,/5). The procedure parallels 
that of Ref. ||^. The resulting RG equations |TB[ follows, 

dyt/dln^= (1 - - Mt)yt + Vji/t 
dyj/dln^= (1 - ej)yj+y1 

dMt/dlni= (2 - - Qyl)Mt (6) 

Here ^ is the running inverse energy scale, yt = t^, and yj = ^J^^. et and ej are the scaling 
dimensions of the t and couplings induced by the on-site interactions V and j|- (where jjj- 
is the longitudinal component of the Jx— coupling), ej < 1 if Jk is antiferromagnetic [§]. et 
can be tuned through 1 depending on the charge-screening interactions @,^]. The parameter 
Mt = F'^Kc^q""''' /2{ac — 1) is generated only when the bare F value is non-zero. The 
scahngs of et, ej, and Ed are not explicitly affected by the Mt term. The mixed- valence 
regime is achieved by tuning the renormalized c?— level Ed to be close to zero. 

Since ej < 1, we can infer from Eq. (^ that yj is a relevant coupling. As for the usual 
Kondo problem, this implies the formation of quasiparticle-like spin excitation spectrum at 
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low energies. The resulting dynamical spin susceptibility has the usual Fermi liquid (^)^ 
form. What determines the nature of the charge excitations in our RG approach is the 
behavior of the yt coupling. If the flow of yt is towards infinity, the charge excitations are 
also quasiparticle-like and a Fermi liquid state emerges at low energies. If the fiow of yt is 
instead towards vanishing or an intermediate coupling, the charge excitations are distinctive 
from the quasiparticle-like behavior; spin-charge separation then takes place in the mixed- 
valence regime. In our case, yj goes to strong coupling logarithmically. This is a sufficiently 
slow fiow towards strong coupling that, when the initial value of tt + Mt is sufficiently larger 
than 1, yt does not fiow towards strong coupling in the regime where the above scaling 
equations are valid. This leads to the possibility for a fixed point with a weak-coupling yt 
despite a strong coupling yj. 

That such a fixed point is indeed a stable one can be seen by a strong coupling analysis. 
We construct a Toulouse-point-like Hamiltonian by first introducing an abelian bosonization 



of the fermionic bath |[T^ and then carrying out a canonical transformation to eliminate the 
longitudinal couplings to these abelian bosons. In the parameter regime for the intermediate 
phase and in the absence of the F— coupling this procedure was carried out in Ref. 0. In 
the presence of the F— coupling, the effective Hamiltonian assumes the following form, 

Hi = ^ E(^^oe~^*^^ + H.c.) - ^(Xn + H.c.) 
V27ra cr 47ra 

H2 = ^(Xn - ^u)(^)(^).=o + i^E(E^- - Xoo){Pk + pl,) (7) 

Here $c,s = ('^'t '^i)/"\/2, where $| and $j are the Tomonaga bosons associated with the 
conduction electron bath, and a denotes an ultraviolet cutoff introduced in the bosonization 
formalism 111. We have taken U = e-^(v^/2)'f^'E.'^^-e^(^/^)*=(S-^---^°°). The Hubbard 



operators are Xqq = |0 >rf ^ < 0|, X^q = \cr >d d < and X^„> = \(t >d d < ct'\fI,F„, 

where and F„ are the "Klein operators" ||T5|-0,|]. Finally, Kg = -\/2[tan~^(7rpo(^ + 
j|^/4)) — tan^-'^(7rpo(^ — is an infinitesimal parameter. The Jj^— term is strongly 

relevant. Its effect is to make the bonding combination of the doubly degenerate | ]>d and 
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I l>d states lower in energy than the anti-bonding combination. The t— couphng leads to 
fluctuations between |0 >d and the bonding combination of the | t>d and | l>d states. It 
renormalizes towards zero since its scaling dimension, 1 + ^F'^Kc^q~°''' , is larger than one for 
any non-zero F. The fixed point with vanishing t and strong coupling is indeed stable. 
Since the t— term is the only coupling that mixes the impurity spin and charge degrees of 
freedom, this establishes the dynamical separation of the charge and spin excitations. 

The strong coupling analysis also allows us to address the effect of the coupling to 
the vector bosons in the effective impurity Hamiltonian Eq. @. The effect of this coupling 
is to add a term 

A^i = ^7E[e^^*^^Ti(0.- + r-^) + H.c] (8) 

k 

where 0^ = (0^ — i(pl)/\/2, to Hi defined in Eq. (|^), and another term 

AH, = gY.{X^^ - XuWk + 0l) (9) 

k 

to H2. In this form, the coupling in AHi is less relevant (in the RG sense) than the 
coupling. On the other hand, AH, leads to an additional screening of the impurity spin by 
the bosonic (p^ bath. 

Consider = 2 first. The scaling dimension of is (/ts/7r)^/2 + {g'^Ks)/2. Given that 
Ks is a vanishingly small parameter, the J^— term continues to have the strong coupling 
behavior for a finite range of g: 

9' <9l = ^ (10) 



The term remains an irrelevant coupling. As a result, the spin-charge separation phe- 
nomenon remains. Physically, while the RKKY interaction introduces mutual screening 
between localized spins and hence inhibits the ability of the conduction electron spins to 
quench the local moments, over this range it is not strong enough to prevent the formation 
of the Kondo singlets. 

While our conclusion that spin-charge separation persists for as = 2 and small g is firm, 
if we extend our strong coupling analysis to large g and assume that AHi remains less 



relevant compared to the J^— coupling even for large g we will predict a phase transition 
for Us = 2: becomes irrelevant for g > g^. This phase transition is entirely the result of 
a competition between the g and Jk coupling. The situation is very different when both g 
and Jk are small. Here we can derive the following RG equations (involving g and Jk only) 
using the standard multiplicative RG procedure |T^], 



dJK/dln^= Jk{PgJk - ^Po-^l: - Ksg"^) 
dg/dlni=-g{]^plJl + K,g^) (11) 



Irrespective of the ratio gj Jk-, Eqs. (|Tl|) imply that Jk always renormalizes towards strong 
coupling. We conclude that, if a phase transition indeed occurs in the spin-isotropic case 
for = 2 it can only happen at a finite g. 



For as <~ 2, we can infer [T^ from Eqs. (0J|,|) a phase transition at g'^^ 



■JkPo 



Ks \ I ■ '^K towards strong coupling for g < gc^a and towards weak cou- 

pling for g > g^^ai both in a power law fashion (in terms of the running scale ^). The 



flows are controlled by an unstable flxed point at Ks/2g* , p^J^'*) = (1, y/2 — as). This 
phase transition has its weak coupling analog. Here, the multiphcative RG equations ( pA]) 
have to be modified by adding a term (2 — as)g to the scaling of g. These RG equations 



then imply a phase transition at gcb = -yr^ '■ Jk fiows towards strong coupling for 
9 < 9c,b and towards zero for g > g^^i,. Again, the fiows are algebraic in terms of ^, and are 
controlled by an unstable fixed point located at {^/K^g* , poJ^) = (2 — a^, 2 — as). Through 
the interference effect described in Eq. (^, an algebraic fiow of Jk towards strong coupling 
will likely make the t— coupling flow towards strong coupling. The latter would destroy the 
spin-charge separation phenomenon as described in this paper. 

One crucial feature is that, in all of the above regimes GcdT) = — < Tt-Co(j{t)cq^{0) > ^^ff 

■imp 

has a long time ^ behavior. This can be seen from -f^e//, Eqs. (^^^, using the bosonization 
representation cq^ = F„^e-^(*=+''*^)/^. 



We are now in position to determine the self-consistent solutions of the lattice Hamil- 
tonian Eq. (|l]). We consider the local interactions such that the solution is the spin- 
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charge separated intermediate phase when the intersite interactions vq = Jq = 0. Here, the 
dynamical spin susceptibihty Xs ~ 72 and the connected dynamical charge susceptibility 
Xc ~ 777 where a is smaller than 2. We ask whether such a solution remains a self-consistent 
one when Jq and vq become finite. Inserting Xs = (e\)^ = (^"ry^ into the self- 

consistency equation (|), the corresponding impurity problem, Eq. (|), has as = 2, Uc = a, 
g = ME^jEs), and F = v^{E^J E,Y'^ (where E^ = [lim^^o - and 
= [lim(^_^o Eg <^(t^ — Wq) / L0°'~^]~^^"'). Also from the self-consistency, GcciT) ~ ^ at long 
times implies that po = Efc S{Ep — Ek) is finite (as already noted earlier). Our analysis above 
then shows that, for Jq smaller than a threshold value {Eg/ E^)gc, the resulting dynamical 
spin susceptibility is also Xs ^ Equally important, the charge susceptibility has the form 
1^ Xc ~ with the same exponent as that of the input Xc- We conclude that, for arbi- 



trary values of Vq and a finite range of Jq, the solution to the lattice model has spin-charge 
separation, with a local spin susceptibility of the Fermi liquid form and a connected 
local charge susceptibility that depends on 1/r through an interaction-dependent exponent. 
The separated spin and charge excitations is most clearly seen in the strong coupling Hamil- 
tonian H^ff given in Eqs. (0, From H^jf we can also infer that the single particle 

Green's function, Gddir) = - < Trdo^{r)dl„{0) >^e//= - < T^(i^,e//(i")4ef f (0) >H,ff, 

imp ^ J J J J 

where c?t,e// — (— l/v^)-^o-o-^CTe~*'^*''/^e*'^*^/^, has a power-law dependence on ^ with an 
interaction-dependent exponent. 

The main remaining issue is the self-consistent solution when Jq becomes sufficiently 
large. What is needed is a method capable of determining from the self-consistent equations 
the evolution of Es towards zero. It is reasonable to assume [^] that as Eg becomes zero 



the local spin susceptibility has the form with < 2. As discussed earlier, in this case 
the spin-charge separation physics will not be likely to hold. Within the spin sector alone 
the situation bears formal similarities with the nontrivial solutions to certain infinite-range 
quantum spin glass problems [pl|-p3[|. 



Our results provide the first step towards establishing spin-charge separation in the two- 
band extended Hubbard model in finite dimensions. The logical next step is to study the 



effects of intersite couplings which are not two-particle in nature. However, in many cases 
the dominant effective intersite interactions are the two-particle RKKY-like interactions we 
have considered. This makes it plausible that the spin-charge separation physics we have 
discussed is already relevant to real materials. 

The weak coupling scaling equations similar to our Eqs. ( ]TT| ) have been independently 
derived by Anirvan Sengupta, who has in addition determined the correlation functions at 
the unstable fixed point for as <~ 2. We would like to thank him for informing us of his 
unpublished results while the present paper was being written. We would also like to thank 
A. Georges, S. Sachdev, and A. Sengupta for useful discussions. Q. S. acknowledges the 
support of an A. P. Sloan Fellowship. 
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